Parametric resonance at the critical temperature in high energy heavy ion collisions 
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Parametric resonance in soft modes at the critical temperature (Tc) in high energy heavy ion 
collisions is studied in the case when the temperature (T) of the system is almost constant for a long 
time. By dividing the fields into three parts, zero mode (condensate), soft modes and hard modes 
and assuming that the hard modes are in thermal equilibrium, we derive the equation of motion for 
soft modes at T = Tc- Enhanced modes are extracted by comparing with the Mathieu equation for 
the condensate oscillating along the sigma axis at T = Tc. It is found that the soft mode of tt fields 
' ■ at about 174 MeV is enhanced. 
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^ . I. INTRODUCTION 

<: 

■ It is expected that a new phase of matter will be formed in high energy heavy ion collisions at RHIC and LHC. 
It is called 'quark-gluon plasma' (QGP) in which the chiral symmetry restoration occurs. Many signals of the chiral 
symmetry restoration have been proposed, but there is no decisive one. One of the proposed signals is the disoriented 
' chiral condensate (DCC) which is a misalignment phenomena in the chiral space. The time development of the chiral 
^ I condensate has been studied in terms of DCC and the possibility of soft mode enhancement by parametric resonance 
' was suggested . Parametric resonance in the last stage of high energy heavy ion collisions was discussed and 
^ , amplified modes were extracted. A parametric resonance is expected even in the chiral phase transition because the 
\^ oscillation of the chiral condensate (momentum fc = 0) may amplify nonzero (k ^ 0) modes. The motion of the 
I condensate must be investigated in order to reveal this phenomena. As indicated in the motion of the chiral 
T-H ' condensate in high energy heavy ion collisions is expected to be almost along the sigma axis in the linear a model. 

The oscillation of the condensate is just that of the sigma condensate. There may be a chance of the amplification of 
the fields by parametric resonance if the sigma field oscillates for a long time. 
O i' The study of parametric resonance at finite temperature in chiral phase transition has never been examined, while 
' that at zero temperature has been performed. Parametric resonance may occur even when the temperature is 
not zero if the condensates moves periodically like a sine function. The motion of the condensates can be described 
(-H ' by the equation of motion with the effective potential if the system is not far from thermal equilibrium. The effective 
potential depends on the temperature of the system which is a function of (proper) time. If the temperature is 
constant for a long time, the condensate will oscillate as it oscillates at zero temperature. The periodicity of the 
motion of the condensate depends on the effective potential at finite temperature. It requires that the temperature 
' must be almost constant in the period of the oscillation of the condensate. 
5^ , It is expected that the temperature may be almost constant as a function of time at the critical temperature {Tc) 
because of the large difference of entropy density between QGP and hadron phases. Especially, it may take some 10 
fm for the phase transition to finish if the phase transition from quarks and gluons to hadrons is of the first order . 
The similar time dependence of the temperature may also occur in the second order chiral phase transition of a high 
energy heavy ion collision. The parametric resonance by the oscillation of the sigma condensate will occur at/near 
Tc in such a case. The similar time dependence of the temperature is also expected when the cooling of the system 
is slow. If the period of one dimensional scaling ]To| is long enough, the temperature is only slowly decreasing as a 
function of time in the last stage of the expansion. 

In this paper, we discuss the possibility of parametric resonance at/near Tc in the framework of the linear a model 
and extract the amplified modes assuming that the temperature is constant at Tc for a long time. The paper is 
organized as follo ws. In Sec. |l^, the equation of motion is derived in the case when the temperature is constant at 
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T = Tc- In Sec. Ill, the amplified modes and the time scale of the amplification are extracted by comparing with 



the Mathieu equation. The time scale is explicitly shown by solving the Mathieu equation numerically. Section IV is 
assigned for conclusions. 
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II. EQUATION OF MOTION FOR SOFT MODES 



The linear a model is a useful tool to describe the motion of the condensate and soft modes below/near the critical 
temperature. The Lagrangian is 



L^\d^W^-\{<^^ -v^f -rH<i,^, (1) 



3 

where (/) — {<pQ, 4'i,4'2,4''i) ~ (f i (/"^ — ^-^d H is the explicit symmetry breaking term. The field (f> is divided 

into three parts, zero mode {(j)jc , the condensation ), soft modes {<j)js) and hard modes {(f>jh)- 

4>j ^ 4']c + 4>]s + 4>jh- (2) 

The bracket, (O), is used to describe O averaged over the hard modes. We apply the free particle approximation for 
the hard modes. Consequently, {4>jh) and are zero if the distribution of hard modes is thermal. Substituting 

eq.(|^) into eq.(|l|) and taking the thermal average of hard modes, we obtain the following effective Lagrangian: 

{C)^{Ck) + {Cv), (3a) 

{Ck) = \d^<Pcd^<f>c + ^d^^sd^'Ps + a^^c^''^. + \{d^'Phd^'^h). (3b) 

{Cv) - -\ {4>l + + 2<l>c ■ 0. - v^Y + H ((/-oc + M 

-^{<Pl) (0c + 0^ + 20c ■ 0s - v') - jii^l)') - (0,, + 0,,)^ (02^), (3c) 



3=0 
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where dot implies the inner product defined by (j>c ■ 4>s = 0jc 

3=0 

If (O) terms have no (pc and d<j)c dependence, Euler-Lagrange equation for (pjc obtained from (£) is 

+ H<ti) i<f>3c + 0,.) + 2A(02J (0^.^ + 0^.^) = 0. (4) 

An equation of the same form for c/jjs is obtained if (O) terms have no 0^ and dcfis ■ Here it is assumed that 
{(f>'jh) ^{T) is j-independent. The meaning of this assumption becomes apparent when the concrete expression of 
J^(T) is obtained in the next section. We introduce the effective potential defined by 

y(0, 0o; T) = ^ (0? + -v^f- Hc^oc (5) 

Note that the order of the phase transition described by this potential is second. 

Since we are interested in the amplification of the soft mode with small amplitude, we first consider eq.(^) with 
(j)js = for all j. This is the zeroth order equation of soft modes: 

□0JC + dV/d^jc = 0. (6a) 

The equation for soft modes with small amplitudes is obtained by substituting eq.(|6a[) into eq.(|^): 

□0,, + A (02 + 02 + 20, ■ 0, + 6T{T) - v^) 4>js + A (20, • 0, + 02) 0,, = 0. (6b) 



Eq.(6b) is the equation of motion for 0s with the back ground field, 0,. 

We would like to calculate the critical temperature. In a realistic case, the critical temperature cannot be defined 
exactly since H is not zero. Nevertheless, one can estimate the critical temperature by requiring that one false 
minimum disappears at T ~ T^. 

(02^ + &F{T,) ~ z;2) 00c - H/\ = 0, (7a) 
0„c = O (n=l,2,3), (7b) 
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where 4>jc is the condensation which is the expectation value of the field (pj at the minimum of the effective potential 
at r = Tc [eq.(|7a|)]. Tc can be obtained by the condition that two solutions of eq.(|7a|) are the same: 



We introduce the fluctuation fields d<j)jc defined by 



1 f4H\ 



6 8 V A / 
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Substituting eq.(||) into eq.(|6a|) and using the condition, eq. ([7a|) , we have 



Ignoring 0{S'^) and higher terms, we obtain 



(8) 
(9) 
(10) 

(11) 



where itiq (iJ/0oc) + 2A0q^ and (77/0oc) [n =1,2,3). It is easily found from eq. (pi]) that the condensate 

oscillates around the minimum of the potential at Tc- In the same way, eq.(|6b|) becomes 



□ + (i7/0Oc) + A + 20c • 4's) 4>js + 2A [(j)c ■ [He] + 4>s ■ (<50c)] 4>. 

+A {ScP^f + A (2^c • 0s + 0^ + 2 {5(j)c) ■ 4's) <Pjc = 0. 



(12) 



Since we are interested in the amplification by the oscillation of the condensate, we consider the small 0^ and discard 
0(0^) and O(0s) terms. The equation for (pjs [eq.(|T^)] in such a case is 



□ c 



0, 



where 



aj =^ i?/0oc + 2A 0c • ((50c) + (0jc)^ + 20JC ((50jc) 



=^ 2A [4>jc4>ic + (fljc (<50ic) + 0ic ((50^0)] = (j 7^ i), 



(13) 

(14a) 
(14b) 



and 0{6'^) terms have been ignored. Since 0oc 7^ and 0„c — 0(n = 1, 2, 3) in the linear a model, the coefficients aj 
and Pji have the following relations: 



P12 = /3i3 = /521 = /523 - /531 - (3: 
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0. 



As stated in Introduction, the condensate moves and oscillates almost along the sigma axis. Then, 
n = 1, 2, 3. Consequently all P is zero. In such a case, eq.(ff3t) is diagonalized: 



The solution of eq. (|11|) for j = is 



□ 00s + [ml + 6A0OC ((500c)] 00s = 0, 

□ 0„s + [m'^ + 2X4>oc ((500c)] 0„s = 0. 

(500c = ~B cos {rriot + 0) , 



(15) 
' for 

(16a) 
(16b) 

(17) 



where B and 9 in eq.(17) are determined by the initial condition of (50oc. This solution, eq.(|l7|), is substituted 



dcf 



into eqs.(16a),(16b) and the transformation from t to — (mot + 0)/2 is applied. The equations of motion now 
become: 



[ de ^ ml \ 



cos(20 0os(e;fc) = o, 



2 , 2 



«^cos(2C);>0„s(e;fc) = o, 



(18a) 
(18b) 



where 0os(C; fc) ^^'^ 0ns (C; k) are the Fourier transformation of 0os('?; x) and 0ras(C; 2?), respectively. 
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III. AMPLIFIED MODES 



Eqs.(| 



( 18b ) are just Mathieu equation. To investigate the amplified modes, we define the following quantities: 

^ (k^ + ml), A^ = ^{k^ + ml). (19) 



The amplified modes are obtained from the above coefficients by the help of the knowledge of Mathieu equation. The 
amplified modes of a field for nonzero modes (fc ^ 0) corresponds to Aa- = 9, 16, • ■ • because Aa- > 4 apparently, while 
that of TT field correspond to A^ = 1, 4, 9, 16, • ■ •. {At^ = 1 is not satisfied for some parameters of the linear a model.) 
Then, the masses of cr and tt fields are needed to determine such modes. The condensate (^oc) at Tc is (4i7/A)^/'^ 
which is easily found from the definition of Tc- Therefore, the masses are obtained if A and H are given. The amplified 
modes corresponding to Aa-, — 1, 4, 9, 16 are shown in Table || for A = 20 and i/^/"^ = 119 MeV which generate ruj^ 
= 135 MeV, m„ = 600 MeV and = 92.5 MeV at T = for v = 87.4 MeV. 

The amplification of tt fields is determined by At^ and the factor (2(5 tt) in the presence of cosine in the Mathieu 
equation. This factor is given by 



2Q. 



SX<j>ocB 



m1\ 



(20) 



where B is the amplitude of zero mode introduced in eq. ( pT| ) . Its numerical value for the previous parameters of the 
linear a model is about 0.051 (MeV~^) x B. For example, 2(3 is —1.53 for B = —30 MeV. Since the solution , w(^), 
of the Mathieu equation for w(0) = l,c?w/(i^(0) = has the property z«(^ + tt) = 
characteristic exponent' which has the relation: 



^w{^) where v is called 'the 



cos(j^7r) — wij:). 



(21) 



The time period, vr, in ^ corresponds to about 2.65 fm in t. It is obtained that |e*'^'^| ~ 3.1 for At^ — 1 and 
—2Qt^/A-p! = 1.53. The numerical solutions for B ^ =f30 MeV for w(0) — 1, dw/dt{Q) ~ are shown in Fig.|l]. lo{^ — tt) 
is about —1.71, which corresponds to |e"^'^| ^ 3.1. It is found that the fields tt are strongly amplified in a few (5 or 
more) fm. 

There are several amplified modes in Mathieu equation, eqs.(18a),( [l8b| ). However, modes larger than k\ cannot be 
amplified modes actually because we assume that hard modes are in thermal. Therefore, the cut off ( fc A) be tween soft 



and hard modes is important in order to know whether the amplified modes extracted from eqs.(18a),(18b) belong to 
soft modes or not. In this paper, /ca is determined as follows. Since 3-{T) is a function of ^a, the latter is obtained 
from eq.(^) in which Tc is related to ^a- Namely fcA is fixed if Tc is given. The density operator of j field for hard 
modes [p^(T)] is assumed as follows: 





/ 








/ 


PhiT) = exp 






/Tr 


cxp 


-T 




I 


\k\>kA J 









T-l \ ^ 



k\>kA 



ujja]{k)aj{k) 



(22) 



where 



mj. We find: 



HT)^{<j^%{x))=i:T[(P,ST)<j,]^{x) 



1 

w 



fc|>fcA 



1 

V 



E 

fc|>fcA 



[exp(w,(fc)/r)-l]c.,(fc) 



(23) 



The first term is the vacuum contribution and the second is the thermal one. The vacuum contribution can be 
discarded in the following calculation since it is the well-known infinity which is removed by the redefinition of the 
energy including the contribution from the soft modes. (Subtraction of the vacuum energy) Eq.(||) is rewritten by the 
massless particle approximation and the integration of angle: 



1 



du — — 

feA/T, exp(u) - 1 



I A ; 



2/3 



(24) 



This approximation is valid near because the masses of a and tt mesons become small. These masses are zero 
at T = Tc in the chiral limit. Note that (O) terms have no dependence on d4>o,(l)o,d4>s and 4>s because ^{T) is 
independent of j in the massless particle approximation. 
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Table |l| shows fcA obtained by solving eq.(|2|) for various with A = 20, H^'^ = 119 MeV and v = 87.4 MeV. 
is about 123 MeV in the chiral limit. It has been estimated from Lattice QCD and some other methods. It is between 
140 MeV and 190 MeV in Lattice QCD (l|Jl|l. Since is above 190 MeV except for < 130 MeV, we conclude 
that fc^ ~ 174 MeV is an amplified mode at least when Tc is adjusted above 130 MeV. Other modes are irrelevant 
since these modes belong to hard modes. ( A:^ ~ 440.1 MeV may be a soft mode in this sense.) 

The amplified mode of tt field is not directly related to the observed one since the mass of tt meson at T ^ 
is different from that at T = 0. Then we must know the effect of the difference of two masses. In quantum field 
theory, the relation between these two masses can be described by Bogoliubov transformation, we can estimate the 
effect by evaluating the coefficients of this transformation. The amplified mode is not changed essentially because the 
difference of the masses between at zero temperature and at the critical temperature is small enough for tt fields in 
this model. The peak coming from the parametric resonance at T = T^, will be found near k^^ ~ 174 MeV if the peak 
is not smeared out by scattering, absorption and so on. 



IV. CONCLUSIONS 



The parametric resonance at the critical temperature in high energy heavy ion collisions is studied in the case 
when the temperature of the system is constant at Tc for a long time. We consider the case in which the condensate 
oscillates along the sigma axis at the critical temperature of the second order phase transition in the framework of 
the linear a model. 

The enhancement of soft mode at about 174 MeV in tt field is found at various Tc above 140 MeV. Other modes 
are irrelevant as amplified modes because these modes belong to the hard modes which are assumed to be in thermal 
equilibrium in the present study. On the contrary, there is no enhanced soft mode in the a field. The amplified mode 
(174 MeV) at T = Tc is softer than that (for example, about 265 MeV in one dimensional expansion case |5|) at T = 
in TT field because the difference between sigma and pion masses at T — Tc is smaller than that at T = 0. It takes a 
short time for the soft modes to be amplified. This implies that the amplification of the soft mode by the parametric 
resonance at Tc is possible in real collisions at high energies. The enhancement by the parametric resonance in the 
present study may be observed experimentally if the smearing effects are weak enough. 

It has been pointed out that there may be the parametric resonance at zero temperature in the last stage of heavy 
ion collisions. Then, the pion momentum distribution caused by the parametric resonance at finite temperatures 
may be the initial distribution of the subsequent parametric resonance at zero temperature. If so, some peaks which 
correspond to the enhanced modes may appear in the final momentum distribution. 

The explicit symmetry breaking term {H4>o) plays important roles for the parametric resonance at T = Tc because 
this term makes the condensate 0oc be nonzero and generates masses at T — Tc - The condensates {(j)jc) and the 
masses {mj) at Tc are exactly zero in the chiral limit {H = 0). Note that eqs. (|l6a| ),(161: ) are valid when O and 
O [(pf) terms in eq. (|l^) are negligible. The resonance structure may change in small H cases. 

We have not shown the strength of the peak in the present study. If we would like to obtain it, we must know the 
time interval in which the temperature is (almost) constant and consider the evolution of QGP. The non-linear effects 
and back reaction are ignored also in this investigation. These problems will be answered in the future studies. 
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TABLE I. Amplified modes in a and tt fields for various and A-^. 



TABLE II. The relation between Tc and fcA. 



FIG. 1. The solution u;(i) of the Mathieu equation for tt fields with = 1 with the initial amplitude B — =p30 MeV and 
w(0) = 1, dw/dt{0) = 0. Thick line and dotted line are for B = —30 and 30 MeV respectively. Dashed line is for —1.71 which 
is the value of cos(f^) at ^ = tt. The arrow with the string '2w/mo' shows the time corresponding to f = tt. 
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k„ fMeV) 


h (MeV) 


1 




174.0 


4 





440.1 


9 


521.9 


682.7 


16 


808.5 


920.5 



Table.| 



(MeV) 


kA (MeV) 


89.5 





130 


150.5 


140 


192.3 


150 


235.7 


160 


280.5 


170 


326.8 


180 


379.3 



Table.g 
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Fig. 1 



